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Abstract 

Structure and dynamics at soft-matter interfaces play an important role in nature and technical 
applications. Optical single-molecule investigations are non-invasive and capable to reveal hetero¬ 
geneities at the nanoscale. In this work we develop an autocorrelation function (ACF) approach 
to retrieve tracer diffusion parameters obtained from fluorescence correlation spectroscopy (FCS) 
experiments in thin liquid films at reflecting substrates. This approach then is used to investigate 
structure and dynamics in 100 nm thick 8CB liquid crystal films on silicon wafers with five differ¬ 
ent oxide thicknesses. We find a different extension of the structural reorientation of 8CB at the 
solid-liquid interface for thin and for thick oxide. For the thin oxides, the perylenediimide tracer 
diffusion dynamics in general agrees with the hydrodynamic modeling using no-slip boundary con¬ 
ditions with only a small deviation close to the substrate, while a considerably stronger decrease 
of the interfacial tracer diffusion is found for the thick oxides. 
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I. INTRODUCTION 


Understanding nanoscale dynamics at soft-matter interfaces is of profound interest in life 
sciences as well as for technical applications. For example, nanoscale structure and dynamics 
are crucial for the performance of organic solar cells [IH3] and inkjet printed electronics [1], 
while, for example, protein adsorption is of great interest in life sciences [5], and in particu¬ 
lar plays an important role in the biocompatibility of implants [H]. The application related 
research is accompanied by fundamental research on the influence of structural changes at 
soft-matter interfaces on material properties^ IE]- Nowadays several well established in¬ 
vestigation methods with nanoscale and even atomic resolution exist, for instance, x-ray 
reflectometry [9] and atomic force microscopy mm- However, their application to the 
investigation of dynamics in soft-matter is still delicate, as they usually demand for sophis¬ 
ticated sample preparation and may cause damage to the materials of interest. The fortune 
of optical investigation methods is their non-invasiveness. Thereby, single molecule methods 
are superior to ensemble methods in retrieving heterogeneous structures and dynamics on 
the micro- and nanoscale [laiis]. Fast diffusion processes and other dynamics related to 
fluorescence fluctuations in the range from nanoseconds to seconds can be monitored by flu¬ 
orescence correlation spectroscopy (FCS) [TEHTB] . Conventional FCS is retrieved employing 
a confocal fluorescence microscope mm- Several more advanced microscope techniques 
have been developed during the recent years [niiiH], including two-photon excitation for 
FCS and also a combination of FCS with stimulated emission depletion (STED) [19]. In 
contrast to many other super resolution techniques, STED does not require for statistical 
averaging, enabling also a high temporal resolution for investigation of fast dynamics [20] . 
For example, it has been successfully applied to discriminate between free and hindered 
diffusion of lipids in a living cell, which had not been possible by conventional confocal 
microscopy m However, drawbacks of STED are its need for alignment of the excitation 
and depletion beams, and a dependence of the spatial resolution on relatively high excita¬ 
tion powers [20], which requires very photostable fluorophores and may damage the sample 
structure. 

Conventional confocal microscopy provides a vertical resolution in the range of 0,5 to 
1 /im [221123], it may be further reduced by using more sophisticated methods [18]. For 
example, 3D STED allows to achieve a vertical resolution of 150 nm [23]. Here we demon- 
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strate a reduction of the vertical resolution in conventional confocal microscopy down to 
50-100 nm by employing optical interference at mirror interfaces. The approach does not 
require any alignment of two or even more laser beams, and it is particularly suited to re¬ 
trieve information about the structure and dynamics at soft-matter interfaces. In principle, 
more advanced methods may also be employed together with mirror interfaces, however, 
the implications of the wavelength dependent interference patterns on such methods have 
to be taken into account. Rigneault and Lenne developed modihcations of 3D correlation 
functions for FCS at mirror interfaces |21H2S] • Below we address the problem of diffusion in 
hlms which are considerably thinner than the focal depth of the excitation prohle. Investi¬ 
gation of thin hlms on this scale is of particular interest for understanding the structure and 
physical properties at soft-matter interfaces. We apply the developed correlation function 
to FCS studies of 100 nm thick liquid crystal (LC) hlms on silicon substrates with varied 
silicon oxide thicknesses d. 

Liquid crystal (LC) materials have been of increasing interest within the recent decades, 
on the one hand due to their abundance in living organisms [271 [28], and on the other 
hand due to an increasing held of technical applications, which does not only span liquid 
crystal displays and optoelectronic devices. For example, recently, LC materials have been 
investigated for enhancement of organic solar cell performance [29] . The here employed 4- 
n-octyl-4’-cyanobiphenyl (8CB) is in its smectic-A state at room temperature and allows us 
to compare the results with literature reports on investigations using diherent experimental 
methods [301 |3T]. LC materials are very sensitive not only to interface structures, but also 
to interaction forces with the substrate, inhuencing their anchoring conditions [321 E3]- For 
the here studied 8CB hlms, in particular, an inhuence of the Si02 layer thickness d on him 
structure and dynamics was found [SI |35]. Here, we vary d in hve steps from native (4 nm) 
to 100 nm, which we expect to shed some light onto this inhuence. 

To allow for a general overview of this study, we start with a short outline of the ex¬ 
periment, before going into details about the developed huorescence correlation function in 
section IIIII and the calculation of the vertical huorescence modulation in section IIVI 
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II. SCOPE OF THE EXPERIMENT 


(a) 8CB dimer (b) PDI 



height in 8CB film z (nm) height in 8CB film z (nm) 


FIG. 1. a,b) Chemical Structures, c) Model (not to scale) of sample structure with 8CB film 
thickness L, Si02 layer thickness d; with refractive indices changing from ny at the substrate (region 
i) to n_L at the interface with air (region iii) and an intermediate region (ii) with mean refractive 
index rim- Orientation of 8CB dimers highlighted as blue cigars, d) Scheme illustrating the use of 
the interference pattern in I{z) for probing the samples, no fluorescence signal from PDI tracers 
(red balls) in dark regions. e,f) Diffusion coefficients according to hydrodynamics: e) Dxy{z) and f) 


Dz{z), calculated from (Cl) and (C2), respectively, for mean diffusion (solid), as well as diffusion 
along (dash) and perpendicular (dot) to the LC director. 


In this work we study thin 8CB dims of thickness L = 110 ± 10 nm on silicon substrates 
differing in oxide thickness d. For details on sample preparation and conduction of mea¬ 
surements we refer to appendix As stated above, the particular structure of the 8CB 
him is very sensitive to interface structure and forces [321135] . In bulk 8CB the molecules 
will arrange paired into dimers with antiparallel dipole moments [36], see FIG. [^a) and the 
LC director oriented parallel to the long axis of the 8CB dimers. We modeled the verti¬ 
cal structure within the 8CB dims using three regions didering by the orientation of these 
8CB dimers, as shown in FIG. [^c), with the 8CB dimers highlighted as blue cigars (not to 
scale). As 8CB is birefrigent, the refractive index depends on the particular orientation and 
is noted also in FIG. [^c), ranging from parallel to the LC director ny close to the substrate, 
to perpendicular to the LC director n± at the interface with air, with an intermediate region 
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for which we employ the mean reflective index rim- Further details on the him structure are 
discussed in section nYi on the vertical huorescence modulation. 

For investigating structure and dynamics in these thin 8CB hlms, we used highly di¬ 
luted perylenediimide (PDI) molecules as huorescence tracers, for the chemical structure 
see FIG. Hf). The PDI tracers are depicted as red balls in FIG. [^), which illustrates 
the huorescence detection depending on the vertical huorescence intensity modulation I{z). 
The dark regions correspond to low values of a typical I{z) caused by negative interfer¬ 
ence. Fluorescence from dihusing PDI tracers can only be detected in the bright regions, 
which depend on d and diher for the two employed excitation wavelengths Xex = 488 nm 
and Xex = 515 nm. Numerical data for each calculated I{z) were employed with the below 
described huorescence intensity autocorrelation function (AGF) to derive vertical and 
lateral D^y dihusion coefficients. Thereby, the vertical dependence of the dihusion coeffi¬ 
cients was modeled by a potential of strength Vq rapidly decaying with distance z from the 
substrate, as will be discussed in more detail in the next section [ITT] and in appendix [B] giving 
details on the derivation of the AGF. 


Since the detection of dihusing tracer molecules is restricted to the areas with sufficiently 
high I{z) (bright areas in FIG. [^), the experimental FGS data report only on the corre¬ 
sponding vertical regions of the him. By incorporating I{z) into the particular AGF for 
htting the data, we take this into account. Gonsequently, the obtained potential and dif¬ 
fusion coefficients are related to the particular vertical region in the hlms spanned by the 
combination of all z for which I{z) > 0.4. As this not necessarily covers the whole him thick¬ 
ness L, we employ a hydrodynamic model with no-slip boundary conditions for calculating 
height dependent dihusion coefficients and comparing them with our results. Details about 
this hydrodynamic model are given in appendix The calculated vertical dependence of 
the dihusion coefficients is plotted in FIG. [^) and f) for Dxy and D^, respectively. Thereby, 
the dihusion anisotropy of 8GB is taken into account by plotting dihusion coefficients along 
(dash) and perpendicular (dot) to the LG director, as well as mean dihusion coefficients 
(solid). 


Before discussing the results from our investigation, we will explain the employed AGF 


and its derivation in the following section III as well as the calculation of / (z) in section IV 
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III. FLUORESCENCE INTENSITY CORRELATION FUNCTION 


In general, for FCS, a temporal (auto-)correlation of the fluorescence intensity (J^) or 
the related (auto-)correlation of the fluorescence intensity fluctuations = J^— < >) 

recorded from a confocal volume in the sample is required [HI IH]. As long as no other 
fluorescence modulation mechanisms in the time region of interest are involved, the diffusion 
coefficients of the ffuorescing probes diluted in the sample may be derived from the ACF. 

The ACF for translational diffusion in FCS experiments is based on the analysis of the 
expression na 


< >=< 5T{t)5T{t') >= J J d^r di^r' I{r) {g{r,t)g{r',t')) I{r'). (1) 

G 


Here, I (r) is the normalized spatial distribution of the detectable light intensity fluctuations 
= /(r)p(r,t)), the function (^(r, t)^(r', t)) is the density-density correlation function of 
the active particles in the detection volume G, and the average is taken over the full time of 
observation. For derivation of the ACF, on the one hand, the particular J(r) in the sample 
has to be known, and on the other hand, knowledge about the expected diffusion process 
is needed. Typically, FCS is used in bulk biological samples, where the spatial behavior of 
/(r) can be approximated by a 3D Gaussian [15], and the density of the active particles is 
expected to be homogeneous. Then the normalized ACF depends on the ratio of the lateral 
Gaussian waist Wxy to the vertical one Wz and has the form [T5| 


^3D(r) 


1 

(1 -F tItd)\/1 + {w^y/WzY{T/TD) 


( 2 ) 


where the correlation time r is scaled with the correlation time for translational diffusion 
te). If FGS is applied to films considerably thinner than the focal depth Wz, the vertical 
diffusion of the ffuorescent probes does not contribute to the AGF. In this case one uses a 
2D Gaussian approximation for the lateral shape of the fluorescence intensity distribution 
/(r) [HEZ], which eventually leads to the 2D AGF 


92d{,t) 


1 

1 + t/td' 


(3) 


In both cases, the AGF (|^, ([^ will yield diffusion coefficients averaged over the focal region 
within the sample during the observation time. 
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By use of a mirror interface, the sampling region can be further reduced due to the 
interference pattern, which for visual light typically yields interference fringes on a scale of 
50 — 100 nm [2H - [25] . Assuming a cosine fringe pattern within the focal depth, a simple 
extension of the 3D ACF (|^ may be used to obtain the correlation function in the form [21] 


9 At) = q^dA) 


1 + A exp- 

D/J 


r 


with Tf = 


A' 




(4) 


where A is a constant, Ag* is the excitation wavelength, n is the refractive index of the 
solvent, and D is the translational diffusion coefficient. Rigneault and Lenne successfully 
applied eq. Q to analyze diffusion of Cy5 in an 80 pm thick water layer on a dielectric 
mirror interface [2S]. They further extended this approach for diffusion in small tubes [2U] . 

Eq. Q can be used as long as homogeneous diffusion within the sample volume is as¬ 
sumed. Structural changes within the sample, in particular close to the interface, as well 
as interactions of the fluorescent probe with the substrate cause heterogeneities EBj- 
If those heterogeneities cannot be neglected, a more precise calculation of the vertical ex¬ 
citation and fluorescence detection is required, together with an advanced model of the 
diffusion of the fluorescent probes. On the other hand, this very interplay between sample 
structure, resulting vertical fluorescence modulation and diffusion coefficients can be uti¬ 
lized to retrieve structure related diffusion coefficients, and thus to increase the amount of 
information obtained from vertical regions in the sample on a scale of 50 to 100 nm. 

Here we propose an ACF approach developed for the diffusion of fluorescent probes in 
thin Elms on mirror interfaces, when the him thicknesses L is considerably smaller than 
the focal depth Wz- Contrary to the previously discussed cases, the vertical huorescence 
intensity distribution I(z) (the zero of z-axis, is placed at the interface of the substrate with 
the LC him, the positive direction is chosen into the thickness L of the LC him) is no longer 


Gaussian, but is modihed by the optical interference (see section IV), causing a vertical 

x^+y^\ point 


‘^'^xy 


modulation of the Gaussians in the x, p-plane, /(r) = /o/(^)exp 
a; = 0, p = 0 is placed at the geometrical center of the laser beam). 

The diffusion takes place in an anisotropic LC which generates a diffusion anisotropy, 
where diffusion parallel to the LC director (Dy) is ~ 1.5 times faster than diffusion per¬ 
pendicular to it {Dj_) [30]. The geometry of the experimental focal volume (see appendix 
0 with 2wxy ~ 0.5 /im > L zz 0.1 pm yields vertical correlation times at least 20 times 
shorter than lateral ones. This allows to distinguish the diffusion processes taking place in 
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vertical and lateral directions. Dne to the diffnsion anisotropy of the material it is feasible to 
employ this experimental possibility and to discern both corresponding diffnsion coefficients, 
the vertical and the lateral diffnsion coefficient D^y The LC director orientation at the 
bonndaries is known from literatnre [39], and one expects that Dz{z) = D\\ and Dxy{z) = Dj, 
as ^ L, while this does not hold for z ^ 0. The interaction of the LC with the snbstrate 
leads to a z-dependence of both Dz{z) and D^y^z), which is modeled by a ^-dependent drift 
term in the diffnsion eqnation, generated by a potential V{z) which rapidly decays with the 
distance from the snbstrate (as z increases), for details abont the mathematical model see 
appendix 

It is clear that the choice of the potential V{z) significantly depends on the geometry 
of the system. According to Blake, diffnsion of a liqnid close to a wall with hydrodynamic 
no-slip bonndary condition can be represented by a potential with a sqnare decay with 
the distance to the wall HOI. For the here stndied thin LC dims, a hydrodynamic no-slip 
bonndary condition at the snbstrate is feasible [TT]. Additionally, the LC properties together 
with the confined geometry give rise to so called psendo casimir forces [12], which can also 
be represented by a 1/z^ decaying potential. Besides the inflnence of the confinement on 
the self-diffnsion of the LC molecnles, the tracer molecnles themselves will interact with the 
snbstrate. For the here nsed nentral perilenediimide tracer molecnles, hydrogen bonding with 
snbstrate snrface silanols is the most likely tracer-snbstrate interaction [38]. According to 
Israelachvili, the potential of hydrogen bonds between two molecnles decays with the sqnare 
of the distance and contains also a dependence on the bond angle [13]. In onr system the 
interface has an nnknown irregnlar distribntion of silanols and thns nnknown bond angles, 
in addition, the interaction potential might be screened by the LC medinm, leading to an 
even faster decay. Long range van der Waals forces acting between the tracer molecules and 
the substrate decay with the cube of the distance [H] . 

As the result of this discussion we model the potential as 

V{z) = V^/z^ (5) 

with some potential strength Vq. Rapid decrease of the potential allows to employ pertur¬ 
bation theory for solving the diffusion equation, for details see appendix Thereby, the 
diffusion coefficients are approximated as being homogeneous within the area of observa¬ 
tion, while their 2 ;-dependence is taken into account by the potential. The other simplifying 


assumptions are that the diffusing isotropic perylenediimide tracers are highly diluted and 
that one can use the reffecting boundary conditions at both interfaces. These assumptions 
are supported by previous single molecule tracking experiments, which in particular showed 
only a small amount of adsorption events (< 17%) [TT] . 

Besides the above explained impacts from the confinement on tracer diffusion at solid- 
liquid interfaces, an influence from the particular LC structure is also expected. In LC 
materials, the interfacial interactions cause distortions of the bulk LC structure which itself 
is induced by the intermolecular interactions. The particular anchoring conditions at the 
substrate depend on the substrate chemistry and roughness [ISHlTj, while the anchoring 
strength also depends on the long range interactions [HI ESI EH]. For thin nematic LC films 
with antagonistic anchoring, i.e. perpendicular (homeotropic) anchoring at the air-interface 
and planar anchoring at the substrate, Lin et ah showed that strong anchoring stabilizes the 
film structure [33], while weak anchoring leads to instabilities of the film thickness and the 
appearance of LC director distortions [HIEU] . For the here studied thin smectic-A LC films 
with antagonistic anchoring conditions, the extent of reorientation at the substrates is still 
not fully understood and thus is not included into the modeled potential (|^. The expected 
influence from the LC reorientation will be discussed along the presentation of our results. 


IV. VERTICAL FLUORESCENCE MODULATION 

To apply the above derived correlation function to the investigation of thin films at 
reffecting substrates, the particular vertical fluorescence intensity modulation has to be 
calculated. Naturally, it is related to the refractive indices of the involved media. Since the 
here employed 8CB is birefringent, its arrangement within the thin film has to be modeled, 
as we briefly stated in section [n| and depicted in FIG. [^c). We now will discuss details of 
the film structure and the employed model before explaining the calculation of the vertical 
fluorescence modulation. 

In the smectic A phase the 8CB dimers form layers spaced by about 3.2 nm aligning 
perpendicular to the layers [36]. At silica substrates 8CB molecules are known to align 
tilted planar, with a tilt angle of ~ 23° to the interface plane 051, while they will align 
perpendicular (homeotropic anchoring) at the interface with air [39] • For smectic liquid 
crystals free energy considerations are strongly unfavorable in respect to bending of the 
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smectic layers [321 [50] . Thus, the arrangement of 8CB in thin hlms with such antagonistic 
anchoring conditions is not straightforward [M]- Just recently, investigations of thin 8CB 
hlms on rubbed PVA were reported [21] • On rubbed PVA, 8CB aligns unidirectional planar, 
in contrast to the randomly oriented planar anchoring at amorphous silica. But also there 
the mismatch of anchoring conditions on the substrate and on the air-interface, causes 
deformations of the smectic LC him. The authors suggest a transition zone to form close 
to the substrate that ranges between 27 and 49 nm into the him, while it may contain a 
number of dislocations or a melted nematic area [2I|. On silicon with native oxide, for very 
thin 8CB hlms a trilayer conformation is reported [2H|- For thicker hlms, the formation of 
a nematic layer at the interface with the substrate is suggested [21 ED]. Recently Roscioni 
et al. performed atomistic molecular dynamics simulations of the 8CB homolog 4-n-pentyl- 
4’-cyanobiphenyl (5CB) at crystal as well as amorphous silica interfaces |17|. For nematic 
5CB on the amorphous silica reorientation from random planar to perpendicular orientation 
occurred within 5 to 15 nm from the substrate in 20 to 25 nm thick hlms m- For our 
considerably thicker smectic LC hlms with L = 110 nm, we expect the reorientation to 
extend further into the him. Preliminary studies of 220 nm thick 8CB hlms on silicon with 
native or 100 nm thick oxide showed an inhuence of the oxide thickness on the him structure, 
suggesting stronger anchoring for 8CB on native oxide [211 [25] and thus a varying extent of 
the reorientations depending on the oxide thickness d. For methodical reasons, we choose 
to use three diherent extensions of the interfacial order spaced by 10 nm, namely 15, 25 
and 35 nm from the substrate into the him to calculate possible huorescence intensities I{z) 
and compare the related ht results. To approximate the random planar arrangement at the 
substrate and the following reorientation, we split the modeled interfacial region into two 
parts: (i) close to the substrate a 3-5 nm thick region with molecules aligned planar and 
refractive index ny, (ii) followed by a region of rearrangement for which we use the mean 
refractive index rim. The remaining top region (iii) of the him consists of smectic layers 
parallel to the substrate with refractive index n±, see FIGj^c). In the following we refer to 
the combination of regions (i) and (ii) as disordered region with height (, motivated by the 
circumstance that at 22° C bulk 8CB is expected to be in its smectic A phase, in contrast 
to this interfacial rearrangement. 

To obtain I{z) we used Essential Macleod"^^ (V.9.7 Thin Film Center, Tucson Arizona) 
to calculate the intensity prohles Ij{(p,^) within our thin hlms for relevant wavelengths j, 
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whereby ^ is the vertical position in respect to the reflecting silicon interface and (p is the 
incident and reflected angle at the film surface. Essential Macleod^^ allows to use uniaxially 
birefringent media within the calculation of reflections in thin Aims. The fluorescent probe 
PDI does not align with the LC, but is isotropically distributed [H]. Within 8 CB its 
fluorescence spectrum contains a peak at Xpeak = 620 nm and a 60% intense shoulder peak at 
Xshoulder 647 Um HU. We model the intensity in emission lemi'-P, as a linear combination 
of the intensity profiles /peak(</3,0> Shoulder0 for Keak, Xshoulder, respectively, and use it 
together with the intensity profiles in excitation lex{y^,0 fo obtain angular fluorescence 
profiles for the two excitation wavelengths, 465 nm and 514 nm according to 


-^( 7^5 0 I^x(SP, 0 [-^peak(7^5 0 T 0-6 fohoulder(7^5 0] 


( 6 ) 


The numerical aperture of 0.9 corresponds to a maximum angle of 64°. Thus, I{z) was ob¬ 
tained by setting z = ^ —d, where d is the oxide layer thickness of the substrate and averaging 
the intensity over the area illuminated by the objective using the approximation 

I{z) = I sin^(0.5°)J(0°, z) + [sin ^(93 -|- 0.5°) — sin ^(93 — 0.5°)] z) 


sin^(64.5° 


ip=i 


(7) 

The above expression is a discrete integral (discretization step ip = 1°) of I{ip,z) over 
the surface of the LC air interface in polar coordinates. The fluorescence intensity from 
a tracer molecule reaching the LC air interface will be decreased by an increasing path 
length of incident and emitted light in the LC medium upon increasing the angle ip 
We approximated this drop of I{p>,z) to be proportional to cos^ 9 ?. Essential Macleod™ 
provides I{ip,^) for s, p and a mean polarization. Here we used the mean polarization, 
because our excitation was circularly polarized, and the isotropic dye molecules could freely 
rotate m The obtained vertical fluorescence intensity modulations I{z) do only contain 
averages over illumination and detection angels. The lateral gaussian intensity distribution 


was already taken into account for the derivation of the correlation function (B2, B4) and 
thus is not included here. 

In close proximity to silicon < 20 nm), I{z) is not only determined by interference, 
but also by non-radiative fluorescence quenching [52]. Several mechanisms for this quench¬ 
ing have been discussed in literature. According to Danos et ah, an exponential distance 
dependence of the ratio of non-radiative to radiative de-excitation will fit most experimental 
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data well [53]. For d < 25 nm we further applied this approximation to I{z) from eq.Q to 
obtain the hnal I{z). 


V. RESULTS AND DISCUSSION 


As outlined above, diffusion dynamics in liquids close to interfaces with solids deviate 
from the bulk behavior. Moreover, for the here studied 8CB dims, the diffusion dynamics 
are related to the - yet not fully understood - local LC structure, and we expect to see an 
influence from the local structure in our results. In this study we compare FCS experiments 
from 10 different experimental conditions, where the optical interference leads to detection 
of the FCS signal from specihc vertical regions in the dims, which are varied by changing 
between two excitation wavelengths Xex and dve oxide layer thicknesses d. This section is 
therefore structured as follows: In part | V A] we demonstrate the induence from the redecting 


substrates on the FCS experiment by comparing the obtained curves from such experiments 
to those on nonredecting glass substrates, and we give a comparative overview of the cal¬ 
culated vertical duorescence intensity modulations I{z). For a detailed discussion of the 
experimental FCS curves for the various combinations of and d as well as the dt results 
from application of the derived correlation function (B2 B4) we refer to appendix]^ In two 
further parts we summarize and discuss the obtained results for the LC structure in the thin 


8CB dims (part VB), and for the didusion dynamics (part VC). 


A. Experimental FCS curves and calculated vertical fluorescence modulations 

A comparison of FCS curves obtained from 8CB dims on silicon wafers to such on glass 
slides shows a faster component (shorter correlation time) for the FCS curves on the sil¬ 
icon wafers, in particular for the native oxide (■), see FIG. b). Diffusion coefdcients 
of perylenediimides in 8CB have been reported to be in the range of 2 to 5 pm^/s from 
single molecule tracking experiments and duorescence recovery experiments at room tem¬ 
perature Neglecting the induence from the interference pattern, a conventional two- 
component dt of FCS curves obtained on silicon wafers yields a one to two orders of magni¬ 
tude larger fast component [30]. Taking the interference pattern into account, the application 


of (B2, B4) yields didusion coefdcients similar to those reported from other experiments. 
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FIG. 2. a) Scaled model (see FIG. visualizing the vertical LG film structure used in c,d). b) 
FCS curves obtained with Aex = 465 nm for PDI in 110 nm thick 8CB films on glass (o) and on 
silicon wafers with native (■) and with 100 nm thick thermal oxide (a). c,d) Galculated fluorescence 
intensities I{z) in 8CB films as a function of distance z to the substrate, using Essential Macleod"*"^ 
for native (blue, solid), 10 nm (purple, dash), 25 nm (pink, dot), 70 nm (yellow, dash-dot), and 
100 nm (red, dot-dash-dot) oxide with c) \ex = 465 nm and d) X^x = 515 nm. 

as will be detailed below. For FCS curves on glass substrates, the small reflections on the 
interfaces can be neglected. However, as can be seen in FIG. [^b), the amplitudes of the 
FCS curves on glass (o) are very low. Conventional FCS is applicable to detect fast diffusion 
down to the range of 1 /im^ /s pSlDE]. Very stable experimental conditions are needed to 
monitor slowly diffusing molecules passing through the focal volume. At the required long 
timescales, background fluctuations compete with the signal form the diffusing dye molecule. 
Consequently, the FCS curve from the 8CB film on glass is on the edge of detectability. 

The situation is different for the films on silicon wafers. There, the vertical fluorescence 
modulation leads to fluorescence fluctuations on shorter timescales. However, the FCS 
curves for the films on 100 nm oxide (a) and on native oxide (■) differ significantly. As 
can be seen from FIG. c,d) the interference patterns highlight different regions in the 
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8CB films on these substrates. For better comparison, in FIG. only I{z) calculated for 
the most reasonable extensions ( of the disordered regions are shown, for more details see 
appendix The dims on native (blue, solid line) and on 10 nm thick oxide (purple, dash) 
were modeled with 3 nm, 12 nm, and 95 ± 10 nm for regions (i), (ii) and (hi), respectively, 
while the him on 25 nm oxide (pink, dot) contained slightly more disorder having extensions 
of 3 nm, 22 nm, and 85 ± 10 nm for regions (i), (ii) and (iii), respectively. For the films on 
70 nm (yellow, dash-dot) and 100 nm (red, dot-dash-dot) oxide, extensions of 5 nm, 30 nm, 
and 75 ± 10 nm, respectively, were used. For a detailed discussion of experimental and ht 
results from the varied combinations of d and X^x we refer to appendix In the following 
part we now will discuss the implications of varying d on the LC structure at the solid-liquid 
interface as found from these results. 


B. Extension of LC reorientation at the solid-liquid interface 

In agreement with previous results [3H |35] we found a structural difference for the dims 
on thin and on thick oxide. Our results point to a varying thickness ( of the interfacial region 
at the substrate within which the LC molecules reorient from random planar |15] to smectic 
layering with the LC director normal to the substrate [M]- For thin oxide {d < 10) nm, we 
hnd C ~ 15 nm, while for d = 70 nm ^ 35 nm. 

The similar surface geometry and roughness for thick and thin oxide [16] induce a sim¬ 
ilar anchoring geometry of the 8CB at the different substrates. But a reduced mobility at 
the solid-liquid interface and a higher anchoring strength are expected for small d due to 
the stronger inhuence from the underlying silicon in that case. Long range forces inducing 
huctuations of the electromagnetic held (as van der Waals forces) are stronger for the more 
polar silicon substrate than for the less polar silicon oxide |l3lllil[5l|. Theoretical consid¬ 
erations [55] as well as experiments on protein adsorption [5l] and dynamics in ultrathin 
liquid hlms [IB] show a decrease in mobility at the solid-liquid interface for an increase of 
long range van der Waals interactions. In LC materials, due to their optical anisotropy, the 
LC structure is related to an anisotropy of the van der Waals interactions [IH] . Additionally, 
the so called ’’Pseudo Casimir” ehect, i. e. thermal huctuations of the LC director held, has 
to be taken into account [12]. All these interactions will contribute to the anchoring strength 
of the LC material at the substrate. Because of the LC reorientation in the him induced by 
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the antagonistic bonndary conditions, a theoretical calcnlation of the anchoring strength is 
not straightforward for the here stndied smectic-A 8CB dims on silicon snbstrates. So far, 
theoretical calcnlations of the van der Waals interaction and the Psendo Casiniir effect, have 
only been reported for simple geometries with no reorientation at the walls nasH]. For 
nematic LC materials, strong anchoring will stabilize the him, while weak anchoring may 
lead to the appearance of structnral changes and instabilities laa ESI SSI Eg. For smectic 
LC dims with antagonistic bonndary conditions, melting of the smectic layers close to the 
snbstrate into the nematic phase is snggested EH ES]. Onr dndings snggest that for thin 
d snch melting exists within a thin region ^ 15 nm, while the remaining part of the dim 

consists of smectic A layers. In contrast, for thick d we fonnd a large extension ( ^ 35 nm, 
of the interfacial disorder in the smectic A 8CB, which matches expectations of an increased 
instability and strnctnral change dne to weaker anchoring fonnd for nematics. For 200 nm 
thick 8CB dims with d = 100 nm, temperatnre cycling into the nematic phase was reported 
to lead to the appearance of strnctnral reorientations, so called focal conic domains (FCD), 
while for dims on native oxide no FCD were fonnd [M]. Althongh FCD were not fonnd on 
spin cast dims withont annealing, their appearance points to a diderent organization of the 
8CB dims on thick and thin oxide. Onr dndings snggest the existence of strnctnral disorder 
close to the snbstrate already for as prepared spin cast dims on thick oxide. 


C. Potentials and diffusion coefficients 

Having discnssed the LC strnctnre within the 8CB dims, we now tnrn to the diffnsion 
dynamics within these dims. An overview of diffnsion coefficients D^y, Dz and potentials 
Vo derived for the particnlar (, is given in TABLE [I] together with the calcnlated vertical 
area of observation for each experimental condgnration. As stated above, the potentials 
V{z) are modeled to acconnt for vertical modnlation of the diffnsion coefficients in the area 
of observation, which not necessarily covers the total dim thickness L. FIG. illnstrates 
Dxy{z) (red) and Dz{z) (bine) using the results summarized in TABLEhowever, for better 
visibility we did not plot the determined errors. Along with the vertical dependence of the 
diffusion coefficients, also the contributing vertical dim regions are illustrated by depicting 
approximations of I{z) as dark and bright regions. We will now summarize and discuss 
the dndings concerning Dxy{z) and Dz{z) from dtting the experimental data. For a more 
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TABLE I. Diffusion coefficients Dxy,Dz and potentials Vb derived from fits to experimental FCS 
curves from 110 nm thick 8CB films on substrates with different oxide thickness d, with modeled 
height C of disordered 8CB at the solid-liquid interface, and indication of vertical region 2 : for 
which I(z) > 0.5 and which thus contributes to the derived diffusion coefficients from htting. 
Errors denote standard deviations from multiple experimental data. 


d [nm] 

Xex [nm] 

Dxy [pm^/s] 

Dz [pm^/s] 

Fo [meV] 

C [nm] 

region 2 ; [nm] 

4 

465 

1.8 ±0.3 

5.2 ±2.6 

-0.006 ± 0.0004 

15 

5 < 2 < 50, 2 > 100 

4 

515 

3.6 ±0.9 

7.5 ±2.9 

-0.006 ± 0.001 

15 

5 < 2 < 50, 2 > 105 

10 

465* 

2.4 ±0.6 

1.8 ±0.9 

-0.007 ±0.003 

15 

2 < 35, 2 > 90 

10 

515 

2.4 ±0.3 

2.3 ±0.5 

-0.006 ± 0.001 

15 

2 < 45, 2 > 90 

25 

465 

2.1 ±0.5 

6.6 ±3.2 

0 

25 

2 > 70 

25 

515 

1.6 ±0.3 

2.6 ±0.9 

0 

25 

2 > 70 

70 

465 

0.04 ±0.02 

0.3 ±0.1 

-0.007 ±0.0005 

35 

2 < 70 

70 

515 

0.2 ±0.1 

1.2 ±0.7 

-0.004 ±0.001 

35 

2 < 80 

100 

465 

0.03 ±0.03 

0.2 ±0.2 

-0.01 ±0.001 

35 

0 < 2 < L 

100 

514 

0.11 ±0.06 

0.3 ±0.2 

-0.006 ± 0.001 

35 

2 < 40 


* values obtained for film thickness L = 120 nm 


detailed discussion of the data and fits we refer to appendix [Pj 

For d = 25 nm, the fluorescence signal is detected from z > 70 nm, see the depiction 
in FIG. 1^ middle row. Within this region, the 8CB is organized in smectic-A layers and a 
vertical modulation of diffusion coefficients should be negligible, in agreement with Vb = 0 
obtained from htting. The derived diffusion coefficients (see TABLE are within errors 
similar to the anisotropic diffusion coefficients D± = 2.4±0.1 pm^/s and P|| = 3.8±0.6 pm^/s 
reported for PDI in 200 nm thick 8CB Elms from single molecule tracking (SMT) experiments 
at room temperature [301 El]. From SMT of a perylenediimide in a 0.5 pm thick 8CB 
cell, anisotropic diffusion coefficients were reported as D± = 2.8 ± 0.5 jjiw?/s and D\\ = 
4.4 ± 0.5 pm^/s IS71. In the thicker cell with parallel alignment at the walls, interfacial 
effects contribute less to the diffusion coefficients, explaining slightly larger values than 
those obtained from our FCS experiment. Both SMT experiments were conducted without 
temperature control, and a slightly larger D may also be explained by the high temperature 
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FIG. 3. Illustration of obtained diffusion coefficients Dxy{z) (red) and Dz{z) (blue) in /rm^/s from 
each experiment for (left) Xex = 465 nm and (right) X^x = 515 nm, and d decreasing from top 
to bottom. For errors in D{z) see TABLE |lj Note the change in scale for thin and thick oxides. 
Approximations of I{z) are depicted within each 8CB film in greyscale as illustrated in FIG. [^). 


sensitivity of 8CB . 

Diffusion coefficients obtained from the films on native oxide are of similar range as 
those for 25 nm oxide, see TABLE On native oxide I(z) contains a high amplitude in 
a narrow upper part of the film (for ^ > 100 nm) and in the lower part 5 < z < 50 nm 
(depiction in FIG. |^top row), whereby the fluorescence signal from the close proximity to 
the interface is weak, see FIG. a) in the appendix As discussed above, C ~ 15 nm 
for films on native oxide. Thus, contributions from reorientation are small and most of 
the signal is recorded from smectic-A layered regions. Therefore, we again can assign the 
obtained diffusion coefficients to the anisotropic diffusion in 8GB as Dxy = D± and Dz = D\\. 
A small attractive potential Vq = —0.006 meV is obtained from fitting, taking into account 
the decrease of the mobility induced by the solid-liquid interface. 

For the remaining three substrates J(0) > 0,5, see FIG. and there is a considerable 
contribution from the area of reorienting 8GB to the detected fluorescence. As a consequence, 
Dxy and Dz can no longer be assigned to diffusion perpendicular and parallel to the LG 
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director. As explained above, within these three substrates we found the extension ( of 
the reorientation region to be smallest for d = 10 nm. For dims on this kind of substrate, 
the fluorescence signal contains contributions from the lower him region {z < 40 nm) and a 
small upper him region {z > 90 nm), see FIG. [^second row. The obtained small potential 
is similar to that for hlms on native oxide, and the obtained D^y ~ 2.4 pm^/s is similar to 
that for d = 25 nm, while is of similar range and thus smaller than for d = 25 nm. This 
may be explained by diherent ehects of the interfacial LC reorientation on dihusion parallel 
and perpendicular to the substrate. Any deviation from the smectic A layering with its LC 
director in 2 ;-direction will lead to a decrease of additional to the hydrodynamic ehect 
represented by the interface potential Vq. On the other hand, if there are extended areas of 
similar orientation with LC director deviating considerably from the substrate normal, D^y 
may even be enhanced. Since D^y is similar to that for d = 25 nm, we expect such areas to 
be small, in agreement with our modeling. 

The situation is different for d > 70 nm, see FIG. [^two bottom rows. As can be seen 
in TABLE [T| hlms on those substrates yielded considerably smaller diffusion coefficients. 
Thereby, the largest values D^y ~ 0.2 pm^/s and Dz ~ 1.2 pm^/s were obtained for d = 
70 nm and X^x = 514 nm (FIG. [^second row from bottom, right). In that case, /(O) ~ 0.85 
is only slightly larger than J(0) ~ 0.7 for d = 10 nm, compare FIG. a) and a) in 
appendix]^ In both cases, I{z) shows large amplitude in the lower him region, containing 
the reorientation of 8CB. For d = 70 nm the huorescence is also detected from the middle of 
the him, decreasing towards upper him regions, while for d = 10 nm, the amplitude of I{z) is 
high again in the upper him region but low in the middle of the him. The attractive potential 
Vo = —0.004 ± 0.001 meV is only two thirds of the potential obtained for d = 10 nm. This 
agrees with the smaller area of observation, which covers two thirds of the him thickness 
only. As can be seen in FIG. [^e,f), hydrodynamic considerations according to (Cl, C2) 
still yield a decrease of dihusion coefficients at 2 ; = 80 nm, which is most pronounced for 
in ^-direction. This ehect can be taken into account for the smaller Dz = 1.2 ± 0.7 pm^/s 
obtained for d = 70 nm and Agx = 514 nm compared to Dz = 2.3 ± 0.5 /rm^/s for d = 10 nm 
(see TABLE [^. Additionally, the suggested larger extension 35 nm for d = 70 nm 
compared to C ~ 15 nm for d = 10 nm will contribute to the smaller Dz, since Dz is highest 
within smectic A layered regions with LC director perpendicular to the substrate. 

In contrast to Dz being still of similar order of magnitude, Dxy = 0.2 ± 0.1 pm^/s is one 
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order of magnitude smaller than D^y = 2.4 ± 0.3 /xm^/s found for d = 10 nm. As discussed 
above, LC reorientation may enhance D^y if it leads to extended areas of similar orientation 
with LC director deviating considerably from the substrate normal. For d = 10 nm we had 
found Dxy to be similar to bulk D±, whereby we concluded that such areas had to be small 
in that case. The decrease of D^y found for d = 70 nm suggests the absence of areas with 
similar oriented LC director on the experimental length scale in the range of 100 nm, in 
agreement with the proposed increased interfacial disorder in the above section. The PDI 
tracer molecules are known to cause a considerable distortion of the local 8CB structure [HU] , 
which in SMT experiments on 200 nm thick 8CB hlms led to smaller diffusion coefficients 
than expected from the comparison with self-diffusion coefficients obtained by NMR [aniiii]. 
In the disordered region close to the substrate, the rearrangement of the surrounding 8CB 
may be even more pronounced accounting for the observed strong decrease of the lateral 
tracer mobility. 

An even stronger decrease of diffusion coefficients is found for d = 70 nm with X^x = 
465 nm (FIG. [^second row from bottom, left) and d = 100 nm with Xex = 514 nm (FIG. 
bottom row, right), yielding one order of magnitude smaller Dz, while the Dxy are decreased 
by one to two orders of magnitude, see TABLE In that case, /(O) ~ 1, decreasing to 
0.5 within 70 and 40 nm from the substrate, respectively. According to the hydrodynamic 


model (Gl, G2), at 70 and 40 nm distance to the substrate there is only a small decrease of 


Di. Hydrodynamic modeling according to (Gl, G2), thus, cannot account for the observed 
strong mobility decrease. However, the signal is recorded from the disordered region with 
little contribution from upper him regions with smectic A layering. As explained above, the 
extended rearrangement of 8GB around the tracer molecules may cause a strong decrease 
in mobility. It is now seen to slow both lateral and vertical diffusion, with still a stronger 
effect on lateral diffusion. For both hlms, potentials similar to those from hlms on thin oxide 
were obtained, see TABLE [T| although here the huorescence is recorded from the lower him 
region mainly, pointing to a stronger vertical modulation of the dihusion compared to that 
for d < 25 nm. 

As mentioned above, for d = 100 nm with Ag* = 465 nm, I{z) > 0.6 in the entire him, 
causing only weak intensities of huorescence huctuations from vertical dihusion, which in 
turn leads to a small amplitude of the FGS curve (see FIG. |^b) rendering further discussion 
of the obtained dihusion coefficients not feasible. 
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In general, stronger long-range interactions with the substrate are expected to cause a 
stronger decrease of the mobility of the liquid at the interface (THl EH ES]- For the here 
studied smectic A dims we find a stronger decrease of the tracer mobility for the films which 
are subject to weaker long-range interactions with the substrate. However, the observations 
made for simple liquids may not hold in the same way for complex liquids as the here 
studied smectic A 8CB. For the latter, a slowdown of tracer diffusion is reported in respect 
to hydrodynamic considerations based on diffusion coefficients for self-diffusion [30], which is 
explained by a necessary rearrangement of the LC around the diffusing tracers. As mentioned 


in the previous section V B, in 8CB dims on 100 nm thick oxide, FCDs representing structural 


rearrangement were found after annealing, while they were not observed on native oxide 
Our dndings suggest that already in the as prepared dims on thick oxide irregularly sized LC 
domains exist. LC materials are very sensitive to surface structure [311 SZj. It is therefore 
reasonable to assume that the grain boundaries are dxed by some surface structure, which 
could be given by irregularly distributed surface silanols known to be present on silicon oxide 
grown on silicon wafers [SB] • On the other hand, transition of grain boundaries by the tracer 
molecules implies their rearrangement, which would cause a slowdown or even repulsion of 
the tracers. 


The absence of FCDs observed for 8CB on native oxide [SI] , points to a diderent internal 
structure of these dims on thin oxides, in agreement with our observations. For thin oxides, 
the underlying silicon enhances the long-range interactions [IHEilEH], thus leading to even 
stronger anchoring of 8CB at the substrate. FCDs as well as dislocations reported for 
8CB on rubbed PVA contain regions with unfavorable perpendicular orientation of 8CB at 
the substrate [SH [39], which seem to be absent in thin 8CB dims on thin oxide. Instead, 
for native oxide, nematic melting of 8CB at the substrate is suggested [SH EO]- In the 
nematic phase, the self-didusion is enhanced in respect to the smectic phase [59], which can 
counterbalance a slightly reduced mobility due to stronger long-range interactions with the 
substrate. Moreover, the absence of grain boundaries for 8CB on thin oxides facilitates the 
tracer didusion in respect to 8CB dims on thick oxides and thus explains the agreement of 
the obtained didusion coefficients with the hydrodynamic model. 
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VI. CONCLUSION 


The here developed correlation function (B2, B4) could be successfully applied to deter¬ 


mine the tracer diffusion behavior of a perylenediimide in thin 8CB dims on silicon substrates 
with varied oxide thickness. In particular, the fit results allowed us to estimate the extension 
( of the interfacial reorientation of the 8CB close to the substrate from three different mod¬ 
eled extensions. We found ( increasing with oxide thickness, which we assign to a decreasing 
stabilization of the liquid crystal structure by long-range interactions with the underlying 
silicon in accordance with literature [161 |33l ESI HU HU IMl ES] . For thin oxides (d < 10 nm), 
C 15 nm was obtained, while for thick oxides {d > 70 nm), a larger disordered region with 
C 35 nm was found. For 200 nm thick 8CB films different LC structures after annealing 
into the nematic phase are reported in literature for native and for 100 nm thick oxides [34] . 
Our results suggest a different structure for the spin cast films even before annealing. 

Our findings not only show a variation in the structural arrangement of the 8CB at the 
substrate. The larger extend of re-orientation for the thick oxides also leads to a stronger 
decrease of tracer diffusion close to the substrate. Thereby, Dz is reduced about one order 
of magnitude (within 50 nm from the substrate), while D^y is reduced even more, about 
two orders of magnitude. This strong decrease of tracer diffusion can be explained by 
dislocations and grain boundaries [3T] in the disordered region, which slow or even repel 
diffusing tracers by hampering the local rearrangement caused by these tracers [SU] . In the 
upper him regions on the thin oxide, anisotropic diffusion coefficients of similar order than 
reported from SMT experiments [301 HUE?] were found, in particular D± = 2.1 ±0.5 pm^/s 
and Zl|| = 6.6 ± 3.2 pm^/s from investigation of dims on 25 nm thick oxide with excitation 
at 465 nm, where ~ 40 nm of the upper him region contributed to the signal. The very 
narrow height (< 30 nm) of this region in case of excitation at 514 nm led to considerable 
errors from htting going along with dihusion coefficients of half of the expected size, which 
we assign to an inhuence from higher orders of the correlation function. Thus, for very thin 


him regions contributing to (B2, B4), higher orders have to be taken into account. 

The huorescence modulation at rehecting substrates may further be applied to investigate 
thin polymer hlms, in particular to yield additional information about the widely discussed 
impact from conhnement ehects to the glass transition temperature [601 El] • Investigations 
of spectral dihusion [62l |63| or rotational dihusion [HU [65| could be studied for diherent 
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vertical regions within snch films withont the need to change the film thickness. 


Appendix A: Experiment 

We performed PCS on liqnid crystal (LC) films on silicon wafers with native (~ 4 nm [54]h 
and thermally grown 10 nm, 25 nm, 70 nm (all: Center for Microtechnologies, Chem¬ 
nitz), and 100 nm oxide (CrysTec, Berlin), employing two different excitation wavelengths 
(Aex = 465 nm and Agx = 515 nm). The silicon wafers were slightly doped (resistiv¬ 
ity 5 — 20 flcm). Glass snbstrates (Carl Roth, Karlsrnhe) were nsed for farther com¬ 
parison. Liqnid crystal 4-n-octyl-4’-cyanobiphenyl (8CB, SYNTHON Chemicals, Wolfen) 
is in its smectic-A phase at room temperatnre and was used without further purifica¬ 
tion. The 8CB films were prepared as previously described mi by spin coating solu¬ 
tions of 50 mg/ml 8CB in Toulene, resulting in film thicknesses L = 110 ± 10 nm. The 
solutions had been doped with N,N’-di-propyl-l,6,7,12-tetra-(4-heptyl-phenoxy)-perylene- 
3,4,9,10-tetra-carboxdiimide (PDI), a perylenediimide dye synthesized by Melanie Bibrach, 
TU Chemnitz, yielding nanomolar concentration of PDI in the 8CB films. Perylenediimides 
are very photostable dyes suited for single molecule investigation. The here chosen PDI does 
not align with the LC structure, but will rotate freely within the LC film, yielding an isotropic 
fluorescence signal on our experimental time scale with correlations times r > 10 /rs BH. 
Chemical structures of 8CB and PDI are given in FIG. [^a) and b). 

To aid focussing, highly diluted Nile Red doped Latex Beads of 0.2 /im diameter (Molec¬ 
ular Probes, Leiden, Netherlands) were spin coated on the clean substrates from solution 
in water before the 8GB film preparation. Defocussing of the confocal excitation on the 
sample in the range of Az >0.1 /im had lead to considerable drops in the amplitude of the 
recorded fluorescence in combination with apparent slower diffusion, as reported by Petrov 
and Schwille in [22]. By aid of the Latex beads, Az cold be kept in a range below 0.1 /rm. 

FGS data were obtained using a previously described homebuilt confocal microscope [H] 
[56] with an avalanche photodiode (SPGM-AQRH-14, Excelitas), a TGSPG-board (Becker 
& Hickl SPG-630) and a microscope objective with high numerical aperture AA = 0.9. 
Since the self-diffusion coefficient of 8GB is very temperature-sensitive |59|, a water cooling 
device together with a ptlOO sensor was implemented into the sample stage. This allowed 
for stabilization at T = 22 ± 0.3°G with a calibration accuracy of ±1°G. The setup was 
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eventually extended by implementing a second excitation beam at 515 nm employing a 25 
mW cw laser (Cobolt Fandango). Circular excitation at the sample stage was obtained 
using sets of A/2 and A/4 plates in each excitation beam [3S]. The excitation light was 
expanded by a linear beam expander (Thorlabs) to £11 the back aperture of the objective. 
To avoid photobleaching, the excitation power was dimmed to 1 /iW (cw equivalent) for 
both excitation wavelengths. Before recording FCS data, a region of 100 /rm x 100 /rm 
on the sample was scanned in reflected excitation light (by removing the long pass Alter 
in the detection beam). Regions showing structural artifacts from preparation could thus 
be excluded from recording. For FCS recordings also a distance of at least 2 /rm was kept 
to Latex beads to avoid the influence from structural rearrangements of the liquid crystal 
close to the beads [301 EZ]. For each of the ten combinations of excitation wavelength 
and substrate a fresh sample was prepared and 10 to 15 FCS measurements were taken at 
different positions on the sample, with an acquisition time of 15 minutes for each one. All 
experimental FCS data sets were fitted separately. For Aex = 465 nm on the 70 nm oxide, 
several data sets contained a rather low amplitude of the autocorrelation curve, pointing to 
some higher experimental instability during measurements. These data sets were excluded 
from further evaluation, leaving only seven valid data sets. In general, the FCS curves were 
less noisy for excitation at 514 nm, which is due to the stronger absorption of PDI at that 
wavelength. 

The focal width Uxy of the exciting laser beam was determined by an approximation for 
diffraction limited focussing [68] 


2Aea; 
'^xy y 

nnA 


(Al) 


N 


where n is the mean refractive index of 8CB at 22°C and Aat is the numerical aperture of 
the used objective. We obtained Wxy = 0.21 /rm and Wxy = 0.23 pm for Ag* = 465 nm and 
Xex = 514 nm, respectively. 


Appendix B: Derivation of density correlation function 

For calculating the density-density correlator entering eq. ([^ we used a mathematical 
model of diffusion in a film of certain thickness with inhomogeneous properties of the media 
along the transverse z and homogeneous ones along the lateral direction. An analytical so¬ 


lution of the diffusion equation is not feasible. As explained in section III, the 2 ;-dependence 


23 




of the diffusion coefficients is induced by the interface with the substrate and decays rapidly 
with distance to the substrate. Thus, for a sufficiently thick him, the 2 ;-dependence of the 
diffusion coefficients is only weak. We use this thick-hlm assumption for deriving approxi¬ 
mated relations for the ACF stated in eq. Q. In more details the model is formulated as 
follows: 


The volume G in eq. Q is the volume of the him of thickness L, 0 < z < L inhnitely 
extended in x and y directions. 


(ii) The properties of the him are isotropic and uniform in the x, j/-plane. Therefore, 
the dependence of the dihusion coefficients can be reduced to dependence on the z- 
coordinate only. Thus the dihusion tensor has only two components, Dxy{z) and 
D,{z). 


(hi) Since the tracers are highly diluted in the him, the density-density correlator (p(r, t)p{r', t')) 
of the active particles can be replaced by their joint probability density, which, in turn, 
can be represented as a product of the conditional probability, P(r,f|r',T), to hnd 
the particle at time t in the small volume positioned at r if it was in r' at time t' and 
of some initial distribution of particles. Ergodicity and stationarity of the dihusion 
process allows to assume that the correlators in ([^ depend only on the time interval 
T = t — t' and can be described by the Fokker-Planck equation of the form 

dP 1 f d^\ 1 d 

supplemented by the initial condition P(r, 0|r', 0) = 5(r — r'). 


(iv) 


There is no hux on the interface of the him. Thus the boundary conditions have the 
form 


a{z)P 


2 ^ 


D,{z)P 


2=0,L 


0 . 


(v) The stationary probability distribution, i.e. the conditional probability P(r,r|r',0) in 
the limit r —)■ cx), depends only on 2 ; and is assumed to satisfy the Gibbs-distribution 
in the potential V{z) (see discussion in section III), 


P(r, r —)• oo|r', 0) 


Nq ^exp 


ksT 
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dz exp 
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The latter being substituted into (Bl) gives rise to an expression for the drift term 


1 d 
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(vi) Since the potential V{z) is generated by the him interfaces, we can assume that in 
case of a thick him it is weak. The same assumption one can make for the dihusion 
coefficients as well and represent them as Dz{z) = Df'* + 5Dz{z) and Dxy{z) = D^y + 
^Dxy (^) • 


The last assumption allows to solve equation (Bl) perturbatively. By the method of 
separation of variables we derive the following expression for the leading order of g{T): 
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Employing a Fourier-transformation in x and y coordinates and a Laplace transformation 

in r, the hrst order correction can be obtained as well: 

'2 “ 

TT 


= 


-'o 


1 + T)i°V/(4M;2 ) 


1 


n=0 


xyJ 

2 2n{0) 

n 7r Dj.y T 


n,m=0 


^xy I 


^ an 


2L2 

an'^O ^ ^ -f- CLni—k 

k=l 




where 


0nm(C) 


g—_ g— 

77,2 — m2 


The Fourier transform Vnm{a,b) of the potential V{z) is given by 
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The zero order term (B2) and hrst order term (B4) now can be used together with 
numerical data for I{z) to derive lateral and vertical dihusion coefficients D^y and as 
well as the strength Vq of the decaying potential for the tracers in the hlms of a given 
thickness L. 
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Appendix C: Height dependent diffusion coefficients according to hydrodynamic 
model 

As stated above, the derived correlation function will only take vertical dynamics into 
account, which occur within the area of observation determined by the vertical fluorescence 
modulation. We therefore compared the results with expectations using continuum hydrody¬ 
namics. Modeling confined diffusion is not straightforward, and analytical solutions for the 
diffusion coefficients derived from hydrodynamic considerations exist only for few cases [69] . 
In general, the proximity of a solid-liquid interface will affect the self-diffusion behavior of 
a liquid, and thus the diffusion of incorporated tracer molecules. Previous experimental 
reports showed a slowdown of the self-diffusion of 8CB at the walls of silica nanopores |7D1 
as well as of the tracer diffusion of perylenediimide in thin 8CB films on silica [H]. For 
the room temperature nematic homolog 5CB, conservation of the bulk diffusion anisotropy 
at the solid-nematic interface was reported cn. Joly et al. discussed the application of 
models from continuum hydrodynamics to nanohydrodynamics. According to their results, 
for nanohydrodynamics on wetting surfaces, a hydrodynamic model with no slip boundary 
condition may be applied [72]. Lin et al. give approximations for the diffusion parallel to 
the substrate D^yiz) and perpendicular D^^z) to it, derived from a simple hydrodynamic 
model employing hard spheres of diameter a and bulk diffusion coefficient D at a solid wall 
with no slip boundary condition 1731 
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Previous experiments have shown that in the here used liquid crystal 8CB the dye 
molecules lead to a large re-orientation of the surrounding LC material [30], which can be 
modeled by an increased hydrodynamic radius a [SU] SI]. In case of PDI diffusion in 8CB, a 
was found to be approximately three times the molecular radius [H], yielding a = 3.5 nm. 
Due to the diffusion anisotropy parallel and perpendicular to the LC director in 8CB, the 
local diffusion coefficient depends on the particular orientation of the 8CB in the differ¬ 
ent regions of the film, as was found for diffusion of CdSe quantum dots in a nematic LC 
cell [71]. FIG.[^e,f) shows D^yi^z) calculated from (Cl) and Dz{z) calculated from (C2) for 
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three different bulk diffusion coefficients obtained for perylenediimide in 8CB, (i) diffusion 
parallel to the liquid crystal director with bulk Dy = 4.4 /im^/s [301EZ] (dash), (ii) diffusion 
perpendicular to the liquid crystal director with bulk = 2.9 /ini^/s [301 EZ] (dot), and 
(hi) a mean diffusion coefficient D = 3.2 pm^/s [331 EZ] (solid). The here studied thin 8CB 
films contain structural reorientations close to the substrate. In that area we expect to find 
deviations from the model. 


Appendix D: Fits to experimental FCS cnrves 

Here we will start the discussion with 8CB films on 10 and on 70 nm oxide, where I{z) 
highlights the disordered regions. Then we will turn to 25 nm oxide, where the fluorescence is 
detected from the upper film region only. On native oxide, there is an additional contribution 
from the lower film region, while the situation for 100 nm oxide depends on X^x- 

From FIG. c,d) it can be seen that for d = 10 nm (dash) and d = 70 nm (dash- 
dot) I{z) has a considerably high amplitude close to the substrate, extending some tens 
of nanometers into the film. As stated above, within this region a rearrangement of the 
8CB molecules from random planar at the substrate 051 to smectic-A layers with the LC 
director aligned along the interface normal should take place |39]. We therefore expect the 
experimental FCS curves to contain some information about this disordered region and its 
extension Q. While for d = 10 nm, there are two short regions (0 nm < z < 40 nm and 
> 80 nm) with high amplitude of I{z), for d = 70 nm, the amplitude is high in one 
large region (0 nm < z < 80 nm). In the latter case, vertical diffusion will cause slower 
fluorescence fluctuations in respect to the situation on 10 nm oxide, leading to an enhanced 
FCS amplitude for slower correlation times for d = 70 nm, which indeed can be seen in the 
experimental curves for correlation times 10“^ s < r < 1 s, compare FIG. |^b) and[^b). 


a. 8CB films on 10 nm oxide - highlighting the disordered region 

For both excitation wavelengths on 10 nm oxide, I{z) is low in the middle of the films 
(z ~ 60 nm) with larger amplitudes for the lower and upper parts of the films, see FIG. |^a). 
The similar I{z) result in similar FCS curves as can be seen in FIG. |^b). The larger noise 
for Xgx = 465 nm (light solid) is due to the weaker absorption of PDI for that wavelength and 
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FIG. 4. Results for 8CB films on 10 nm Si02. a) I{z) calculated for Xex = 465 nm (blue) 
with C = 15 nm (solid), Q = 2b nm (dash), and C = 35 nm (dot), as well as for X^x = 514 nm 
(green) with C = 15 nm (dash-dot), (" = 25 nm (dot-dash-dot), and (" = 35 nm (dash-dot-dash), 
b) Experimental FCS curves for Xex = 465 nm (light solid) and Xex = 515 nm (black solid) with 
fits for Xex = 515 nm and L = 110 nm and ^ = 15 nm (red, dash), ^ = 25 nm (blue, dot) and 
C = 35 nm (yellow, dot-dash-dot), c) Fitting errors for Xex = 465 nm (♦) and Xex = 515 nm 
(0)- d,e) Diffusion coefficients Dxy (♦) and (□) obtained from fits, c-e) Values obtained for 
L = 100 nm, L = 110 nm, and L = 120 nm displayed with 1 nm spacing. Error bars denote 
standard deviations using multiple data sets. 


results in a larger scatter of Dz, see FIG. |^d, □), because of the shorter correlation times 
related to vertical diffusion. For Xex = 515 nm best hts were obtained for L = 110 nm and 
C = 15 nm (FlG.|4|b, red dash). Fits with (^ = 25 nm (blue dot) and C = 35 nm (yellow dot- 
dash-dot) showed slightly larger deviations. This matches the situation for Xex = 465 nm, 
if only the part of the FGS curve for r > 1 ms is taken into account (hts not shown). 


For the upper him region {z > 80 nm), the him consists of smectic A layers, which leads 
to Dxy = D_!z < D\\ = Dz- However, in the disordered region at the substrate, the dihu- 
sion anisotropy should lead to Dxy > Dz, which may even be enhanced by hydrodynamics 


showing a stronger decrease of the vertical dihusion coefficient, compare (Gl) to (G2). The 


hrst peak of I{z) contains the disordered region, but foYz>C. also smectic layers parallel to 



























the interface, where D^y < D^. As can be seen in FIG. |^a), for increasing ( the amplitude 
I{z) for the hrst peak {z < 40) nm decreases, while the amplitude for the upper him region 
{z > 80) nm stays almost the same. Therefore, the contribution from the disordered region 
to the FCS curve decreases with increasing (, and we would expect D^y ~ for (^ = 15 nm, 
but Dxy < Dz for larger C. For Aea; = 515 nm and for (^ = 15 nm, D^y (♦) and (□) 
are approximately equal, matching the expectation, see FIG. |^e). However, for ( = 25 nm, 
D^y > Dz contrary to the expectation, apart from ^ = 35 nm, where Dz is slightly larger 
than Dxy for L = 120 nm. But for this C, fhe amplitude of I{z) in the disordered region 
is very small and Dz should be considerably larger than Dz;y. Thus, the best hts and the 
obtained diffusion coefficients both point to the size of the disordered region to be in the 
range of 15 nm. The variation of him thickness L within htting had only minor ehects in 
case of Xex = 515 nm, while for X^x = 465 nm most reasonable dihusion coefficients were 
obtained for L = 120 nm pointing to a slightly thicker him for that particular sample. 


b. 8CB films on 70 nm oxide - highlighting the disordered region 

As can be seen in FIG. a), on 70 nm oxide, I{z) > 0.5 for z < 80 nm, with a shift 
of the peak to smaller z for Xex = 465 nm (blue) compared to Xex = 514 nm (green). 
Gonsequently, the FGS curves obtained at Xex = 465 nm contain a larger contribution for 
longer correlation times r, related to slower dihusion close to the substrate, see FIG. [^b). 
Best hts were obtained for (^ = 35 nm, as can be seen in FIG. |^b) showing an experimental 
FGS curve for Xex = 515 nm (black solid) together with hts for L = 110 nm and C = 35 nm 
(yellow, dot-dash-dot), ( = 25 nm (blue, dot), and C = 15 nm (red, dash). ^ = 35 nm also 
yielded the smallest errors with the most narrow standard deviations, see FIG. [^c). 

The high amplitude of I{z) for z < 80 nm covers the disordered region at the solid- 
liquid interface, as well as some part of the adjacent him region organized in smectic A 
layers. For both excitation wavelengths, I{z) calculated for ( = 35 nm yields the highest 
relative contribution from the disordered region, see FIG. [^a) (blue dot) for Xex = 465 nm, 
and (green dash-dot-dash) for Xex = 514 nm. As discussed previously, in the disordered 
region we expect Dz < Dxy and in the smectic A layers Dz > Dxy, which here should lead 
to Dz > Dxy, due to the larger contribution from the smectic A region. Fitting results 
match this expectation, see (FIG.[^d,e), with only minor inhuence from variation of L. For 
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FIG. 5. Results for 8CB films on 70 nm Si02. a) I{z) calculated for X^x = 465 nm (blue) 
with C = 15 nm (solid), C = 25 nm (dash), and C = 35 nm (dot), as well as for X^x = 514 nm 
(green) with C = 15 nm (dash-dot), C = 25 nm (dot-dash-dot), and C = 35 nm (dash-dot-dash), 
b) experimental FCS curves (solids) with fits for Xe,x = 515 nm. c) fitting errors for Xex = 465 nm 
(♦) and Xex = 515 nm (0). d,e) diffusion coefficients Dxy (♦) and (□) obtained from fits, 
c-e)Values obtained for L = 100 nm, L = 110 nm, and L = 120 nm displayed with 1 nm spacing. 
Error bars denote standard deviations using multiple data sets. 


Xex = 514 nm, most reasonable diffnsion coefficients are obtained for (^ = 35 nm, agreeing 
with best fits seen in FIG. [^b) (yellow, dot-dash-dot). For Xex = 465 nm only 7 FCS cnrves 
conld be analyzed, the resnlts follow the trend seen for Xex = 514 nm. 


c. 8CB films on 25 nm oxide - upper film region only 

For d = 25 nm all calcnlated I{z) show a signihcant amplitude only for z > 70 nm, 
see FIG. |^a). Consequently, the fluorescence can be detected only from the upper part of 
the 8CB film leading to similar experimental FCS curves for both excitation wavelengths, 
see FIG. b), and the £t results do not depend on see FIG. c-e). The rather large 
fitting errors are induced by the narrow vertical film region contributing to the signal. On 
short vertical distances, higher orders of the correlation function might have to be taken 
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FIG. 6. Results for 8CB films on 25 nm Si02. a) I{z) calculated for X^x = 465 nm (blue) 
with C = 15 nm (solid), C = 25 nm (dash), and C = 35 nm (dot), as well as for X^x = 514 nm 
(green) with (" = 15 nm (dash-dot), (" = 25 nm (dot-dash-dot), and (" = 35 nm (dash-dot-dash), 
b) experimental FCS curves (solids) with fits for Xex = 515 nm. c) Fitting errors for Aex = 465 nm 
(♦) and Xex = 515 nm (0). d,e) Diffusion coefficients Dxy (♦) and Dz (□) obtained from fits, c-e) 
Values obtained for L = 100 nm, L = 110 nm, and L = 120 nm displayed with 1 nm spacing. 
Error bars denote standard deviations using multiple data sets. 

into account for improving the fits. The only small observed vertical area also leads to large 
standard deviations for the Dj, see FIG. |^d,e), which are most pronounced for X^x = 465 nm 
and Dz^ due to the large noise at short r in that case (FIG. [^b)). Due to the only small 
upper vertical region contributing to the signal, the actual him thickness L has a considerable 
inhuence on the obtained D{z), as can be seen in FIG[^d,e). 


d. 8CB films on native oxide - lower and upper film regions 

Similar to the situation for d = 10 nm, on native oxide I{z) contains two peaks for z < L, 
but there is a considerably smaller contribution to I{z) from the close proximity to the solid- 
liquid interface {z < 5 nm), see FIGj^a) and thus from the disordered region. We therefore 
expect the diffusion coefficients to be larger than for d = 10 nm, with Dz > Dxy Also here 
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FIG. 7. Results for 8CB films on native Si02. a) I{z) calculated for Xex = 465 nm (blue) 
with C = 15 nm (solid), ^ = 25 nm (dash), and C = 35 nm (dot), as well as for Xex = 514 nm 
(green) with C = 15 nm (dash-dot), (" = 25 nm (dot-dash-dot), and (" = 35 nm (dash-dot-dash), 
b) experimental FCS curves (solids) with fits for Xex = 515 nm. c) fitting errors for Xex = 465 nm 
(♦) and Xex = 515 nm (<>)• d,e) diffusion coefficients Dxy (♦) and (□) obtained from fits, c-e) 
Values obtained for L = 100 nm, L = 110 nm, and L = 120 nm shown with 1 nm spacing. Error 
bars give standard deviations using multiple data sets. 


best fits are obtained for (^ = 15 nm, as can be seen exemplarily in FIGj^b). Fitting errors 
are similar within standard deviations from multiple data sets (FIG. [^c). For Xex = 514 nm 
and C = 35 nm we obtained Dxy ^ D^, which is unreasonable. For ( = 25 nm, Dxy ~ D^ 
within standard deviation, while only for (^ = 15 nm, D^y < D^, as expected. Thus, also on 
native oxide, best hts and the obtained diffusion coefficients both point to C = 15 nm. The 
larger standard deviations obtained for D^, see FIG. [^d,e), are caused by the noisy FGS 
curves at short r for both Xex, see FIG. [^b), and the narrow 2 ; > 100 nm where I{z) > 0.5. 


e. 8CB films on 100 nm oxide - highlighting the solid-liquid interface 

As can be seen in FIG. [^a), for d = 100 nm the general shape of / (z) obtained for different 
( is similar for each excitation wavelength, which leads to similar htting errors (FIG. He)) 
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FIG. 8. Results for 8CB films on 100 nm Si02. a) I{z) calculated for X^x = 465 nm (blue) 
with C = 15 nm (solid), C = 25 nm (dash), and C = 35 nm (dot), as well as for X^x = 514 nm 
(green) with C = 15 nm (dash-dot), (" = 25 nm (dot-dash-dot), and (" = 35 nm (dash-dot-dash), 
b) experimental FCS curves (solids) with fits for Xe,x = 515 nm. c) fitting errors for Xex = 465 nm 
(♦) and Xex = 515 nm (0). d,e) diffusion coefficients Dxy (♦) and (□) obtained from fits, note 
the smaller scale, c-e) Values obtained for L = 100 nm, L = 110 nm, and L = 120 nm displayed 
with 1 nm spacing. Error bars denote standard deviations using multiple data sets. 

and only minor variations in diffusion coefficients (FIG. |^d,e), zoomed in scale!), rendering 
a discrimination of C, from this experiment not feasible. 

A comparison of calculated I{z), see FIG. [^c,d), yields the highest influence from the 
solid-liquid interface for d = 100 nm. This is most prominent for Xex = 514 nm, because 
here I{z) decreases to about 40% of its initial value within 40 nm distance to the substrate 
(FIG.ga, dash-dot-dash). In contrast, for Xex = 465 nm and in particular for (^ = 35 nm, 
I{z) > 0.6 in the entire film (FIG. a, blue dot). This leads to an only small intensity 
of fluorescence fluctuations for vertical diffusion, which in turn causes a small absolute 
amplitude of the FGS curves, see FIG[^b, A), approaching the situation on glass substrates 
(FIG 1^ b, o). The smaller influence from vertical diffusion also leads to an increase in 
the relative amplitude for longer r in the normalized FGS curve obtained at Xex = 465 nm, 
compared to the curve at Xex = 514 nm, see FIG.j^b) light solid and black solid, respectively. 
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For \ex = 465 nm, fitting renders D^y and Dz of similar size as their standard deviations 
from multiple fitting, therefore, a further discussion is not feasible. 
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